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Abstract 

It is shown that gravitational theory allows stable equilibrium con- 
figurations of the degenerated Fermi- gas, whose masses M ^ Mq 
and sizes are less than a = 2GM/c ( G is the gravitational constant 
and c is the speed of light). They have no events horizon at the dis- 
tance r = a from the center. This are objects with low luminosity and 
are a new candidate in Dark Matter im the Universe. Orbits of test 
particles near the above objects are considered. 



1 Introduction. 

Our knowledge about the limiting masses of equilibrium stars configu- ra- 
tions and about the fate of collapsing objects is based on the predict- ions of 
Einstein's general relativity. The astrophysical data do not contradict these 
results, however, they do not prove their validity unambiguously. When an- 
alyzing the observation data it is necessary to take into consideration the 
fact that there also exist other gravitation theor- ies which satisfy (like the 
general relativity) the post-Newtonian tests, predicting, however, other phys- 
ical results in strong fields. (See, for example , [Q] and ). In this paper, 
starting from the gravitational equations proposed before , we show that 
stable equilibrium configurations of a degenerated Fermi-gas can exist with 
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the following properties: their masses M are many times as large as the Sun 
and their sizes are smaller than a = 2GM/c , where G is the gravitational 
constant and c is the speed of light. 



2 Spherically-Symmetric Gravitational Field. 

The above mentioned gravitation equations are 
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These are vacuum bimetric equations for the tensor 



(Greek indices run from to 3) , where 
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Tap are the Christoffel symbols of the pseudo-Euclidean space-time E whose 
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fundamental tensor is rjap , F^^ are the Christoffel symbols of the Riemannian 
space-time V of dimenstion n , whose fundamental tensor is gap ■ The comma 
in eqs. (|ID denotes the covariant differentiation in E. 

A peculiarity of eqs.(P is that they are invariant under arbitrary trans- 
formations of the tensor gap retaining invariant the equations of motion of a 
test particle , i.e. geodesies in . In other words, eqs.(|I|) are the geodesic- 
invariant . Thus , the tensor field gap is determined up to geodesic mappings 
of space-time V (in the way analogous to that the potential A in electrody- 
namics is determined up to gauge transformations). Therefore, additional 
conditions can be imposed on the tensor gap ■ In partucular, |^ if the tensor 
gap satisfies the conditions 
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then eqs (|I]) will be reduced to vacuum Einstein equations Rap 
Rap is the Richi tensor. As distiguished from the gap or F^^ the tensor B^^ is 



, where 
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invariant under geodesic mappings of space-time V. It is in a sense analogous 
to strength tensor Fa/B in electrodynamics which is invariant under the gauge 
transformations. Let us find the spherically symmetric solution of eqs.(|^) 
for a distant observer . We suppose that the observer's frame of reference 
is an inertial one and space-time is pseudo-Euclidean. Choose a spherical 
coordinate system. Then, if a test particle lagrangian is invariant under the 
mapping t — > — t , the fundamental metric form of space-time V can be 
written as 

ds^ = -A{drf - B[{def + sm\e){d^f] + C{dxy, (6) 

where A,B and C are the functions of the radial coordinate r . Proceeding 
from the above- stated, we shall find the functions A ,B and C as the solution 
of the equations system 

Ra^ = 0. (7) 

and 

Qa = 0, (8) 

which satisfy the conditions: 



lim ^ = lim 5 = lim B = 0. (9) 

r— >oo r— >oo r— >oo ^ ' 

The equations Ru = and Rqq = can be written as 

BL2 - 2CLi = (10) 



and 
where 



B'C - 2BL2 = 0, (11) 



L, = i?i2i2 = B"/2 - {B'f/m - B'A'/{4A), 

L2 = i?ioio = -C"/2 + C\AC)'/{AAC) 

and the differentiation oi A , B and C with respect to r is denoted by primes. 
Eqs.(ll) give 

B^AC = r\ (12) 
First, the combination of eqs. ([T0|) and ( pT]) gives 

ACLi - B'C = 0, 
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I.e. 



2B" - {B'Y/B - B'{ACy/{AC) = 0. (13) 
Also, taking the logarithm of eq.(0) we obtain 

{AC)'/{AC) = 4/r - 2B'/B = 0. (14) 

Equation (p!3| ) then becomes 

2B" + {B'Y/B - AB'/r = 0, 

or 

u'/u = 4/r, (15) 
where u = {B'^B . By using (H) we find 

B = {r' + f3Y\ (16) 

where /? is a constant. 

Next, from eq. we find by using eqs.(P) 



C = l- a/B'/\ (17) 

where a is a constant. The non-relativistic limit should be considered in 
order to find this constant . The equations of motion of the particle will 
coincide with the Newtonian ones if a = 2GM/c^ . 

Finally, we can find from eq. ([T^) the function A. The constant j3 in 
eq.(p!6|) remains indeterminate. Since Qa is a vector, the gauge condition (H) 
is covariant. Different constants f3 cause the differnt solutions of eq.(|l]) in 
the used coordinate system . Setting f3 = a we shall obtain the functions A 
, B and C which have no singularity at r ^ : 

A = {ff{l - a/f)-\ B = f,C=l- a/f, (18) 

where 

/ = (r3 + a3)i/3 

and f = df/dr . 

We consider solution (|T^) as the basis for the analysis below. 
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3 Equilibrium configurations . 

The gravity force in the spherically - symmetric field is given by 

= -m{Bi^ - Blpx')x''x'^ (19) 

where i;* = dx^/dt .(Latin indices run from 1 to 3). The gravitational 

force affecting the a test particle of the mass m in rest ( = ) is given by 

GmM . , , , 

1 - «//) , (20) 



r 
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where G is the gravitational constant, a = 2GM/c^, c is the speed of light, 
/ = {a^ + r^y^^ . Fig. 1 shows the plot of the function Fi = —(1/2 r^)(l — 
a//) against the distance r = r/a. 
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Fig.l The plot of the function Fi against the r = r/a 



It follows from Fig.l that the \F\ reaches its maximum at r of the order of 
a and tends to zero at r — > 0. It would therefore be interesting to khow what 
the limiting masses of the equilibrium configurations the gravitational force 
F{r) ( pOD can admit. To answer this question we start from the equation 

dp GpM 



dr r 
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(21) 
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In this equation p is the pressure , M = M{r) is the matter mass inside of a 
sphere of the radious r, p = p{r) is the matter density at the distance r from 
the center, a and / is the function of M(r). 

Suppose the equation of state is p = Kp^ , where K and F are constants. 
For numerical estimates we shall use their values : 

For a degenerated electron gas: 

F = 5/3 K = 1 ■ 10^^ SGS units at p ^ po? where po = 10^ gm/cw?, 
r = 4/3 K = 1- 10^5 SGS units at p > po- 
For degenerated neutron gas: 

F = 5/3 K = 'b ■ 10^ SGS units at p -C po, where po = 5 ■ IQ^^gm/cw? , 
r = 4/3 K = 1- 10^5 SGS units at p > po- 

For rough estimates we replace dp/dr by —p/r, where p is the average 
matter pressure and R is its radius. Under the circumstances we obtain from 
eq.lgj) 

^-^(1-//). (22) 

U a , then the term a/ f is negligible. Setting M{R) ^ pR^ we find 
the mass of equilibrium states as a function of p : 

M=iK/Gf''p^^-^/''^^'/^\ (23) 

It follows from eq.(p3|) that there is the maximal mass ^ M = (K/G)^^'^ 
at p > Po . 

However, according to eq.(^) , there are also equilibrium configurations 
at R < a . In partucilar, at i? -C a we find from eq.(^) that the masses of 
the equilibrium configurations are given by 

M = c9/2l0-lir-3/4G^3/2^-(r-l/3){3/4)_ (34) 

These are the configurations with very large masses. For example, the fol- 
lowing equilibrium configurations can be found: 

the nonrelativistic electrons: p = lO^gm/cm^, M = 1.3 ■ lO'^'^gm, R = 
2, 3 ■ lO^^cm, 

the relativistic electrons: p = 10''^ gm / cm^ , M = 2.3 ■ 10'^^ gm, R = 1.3 ■ 
lO^^cm, 

the nonrelativistic neutrons: p = lO^'^gm/cm'^ M = 3.9 ■ lO^^gm, R = 
1.6 ■ lO'^cm. 
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The reason of the two types of configurations existence can be seen from 
Fig. 2, where for p = 10^^ gm/cm? the plots of right-hand and left-hand 
sides of Eq.(3) against the mass M are given. 
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Fig2 The plot of right-hand (1^2 (M)) and left-hand {Wi{M)) sides of 
Eq.(|2|) against M. 

The following conclusions can be made after considering the plots of the 
above kind: 

1. There are no equilibrium configurations whose the density is larger 
than a certain value pmax ~ lO^^gm/cm^ . 

2. For each value of p < pmax there are two equilibrium states (with 
R> a and R < a). 

Are the configurations with largre masses stable? 

The total energy of the degenerate gase is E = Eint + Egr, where E^t is 
the intrinsic energy and Egr is the gravitational energy. The gravitational 
energy of a sphere is 




R 



(25) 



where 



X(r)= Fdr' (rO-^(l-l//), 

J oo 



7 



a = 2GM{t)/c\ f = {a{rf + (r')')^/^ 

M(r) = 47r dr' p ir'f. 
Jo 

The function x(r) is approximately 

X{r) = (l/r)(l - exp(-r/«)). 
Therefore, at p = const up to a constant of the order one 



E = 



R 



-(1 — exp(— i^/a). 



(26) 



(27) 



The intrinsic energy Ei^t = J u dM, where u is the energy per the mass 
unit. For the used equation of state u = K{T — l)~^p^~^ . Thus, up to a 
constants of the order of one 



E = KMp^-^ - GM^'^p^/^[l - exp(-gM-2/3p-V3]^ 



(28) 



where Q — /2G . As an example. Fig. 3 and Fig.4 show the plot of the 

function E = E{p) for the nonrelativistic neutron gas of the mass M — 
10^^ gm and M = 10^^ gm (neutron stares) correspondingly. 
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Fig. 3 The plot of the function E = E(p) for the neutron configuration 
of the mass M = 10^® gm. 
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Fig. 4 The plot of the function E = E{p) for the neutron star of the mass 
M = 10^^ gm. 

The analysis of such plots show that the function E = E{p) has the 
minimum. Thus, the above equilibrium states of large masses are stable. 

4 Orbits of non-zero mass particles. 

The equations of motion of a test particle of non-zero mass in the spherically 
symmetric field are given by 

= {c^C/A)[l - {C/E){1 + a'l^/B)], (29) 

= c CJa/{BE) (30) 

where (r, ip, 6) are the spherical coordinates ( is supposed to be equal to 
7r/2) , r = dr/dt , ip = dip/dt , E = E/{mc) , J = J/{amc) , E is the particle 
energy, J is the angular momentum. Let u = 1// , where / = (1 -\-r^y/^ 
and f = r/a . Then the differential equation of the orbits , following from 
eqs. ([29| ) and , (pOD can be written as 

du/d(p = Q{u), (31) 
9 




where 

Eq.(^) differs from tlie orbit equations of general relativity by tlie 
function u instead of the function 1/r . Therefore, the distinction in the 
orbits becomes apparent only at the distances r of the order of 1 or less than 
that. 

Setting r = in eq. ( ]29|) we obtain E = Af{r) , where 

M{r) = {l-l/J){l + j'/f) (32) 

is the effective potential . Fig. 5 shows the function A/" = A/'(r) at the 
typical value of J = 4 . 
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Fig. 5 The effective potential Af = Af{r) at J = 4. The horizontals E = Ei 
and E = E2 show the types of particle orbits due to lack of the event horizon. 



The function 7V(r) differs from the one in general relativity in two re- 
spects: 1) it is determined at every point of the interval (0, 00) and 2)it 
tends to zero when r ^ 0. 

Possible orbit types can be shown by the horizontals E = Const . The 
orbits with energies Ei and E2 in Fig. 5 have the pecularity due to the lack 
of the event horizon. The first orbit type ,with energy Ei , exceeding the 
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maximum of the function z/(r) , shows that for each value of J there exists 
such value of E that the gravitational field cannot keep a particle escaping 
from the center. Fig. 6 shows the difference in the orbits of the above type 
and the ones in general relativity. 
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Fig. 6 The particle orbits ak E = 1 and J = 1.99 accordig to eq.(0) (curve 
1) and the one in general relativity (curve 2, dots) . At the distancees r > 1 
the curves are very close together. The coordinates x and y are determined 
in the following way : x = fcos{ip) and y = rsm{ip). 

The orbits of the second type, with the energy E2 , are the ones of particls 
kept by the gravitational field near the field center. 

It follows from eqs. (p9|) and (|30|) that the velocity of a test particle freely 
falling to the point mass M tends to zero when r ^ . The time of the 
motion of the particle from some distance r = tq to r = is infinitely large. 

In the effective potential plot the stable circular orbits are shown by the 
minimum points of the function A/'(r) , the unstable ones — by the maximum 
points .The minimums of the function J\f{r) exists only at J > which 
corresponds to the value of the function /(r) > 3 . Therefore, stable circular 
orbirs exist only at r > r^r , where r^r = ce ~ 2.96 a . The orbital 

speed of the particle with r = r^r is equal to 0.4c . At r < unstable 
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circular orbits can exist. At J ^ oo the location of the maximums tends 
to / = 3/2 . Therefore, the minimum radius of the nonstable circular orbit 
is fmin = 1.33q! . (In general relativity it is equal to 1.5a) .The speed of 
the motion of a particle on this orbit is equal to 0.51c. The binding energy 
E = 0.0572 , just as it occurs in general relativity. 

The rotation frequency uj = cp of the circular motion will be 

uj=[c(f~l){Ta)](J/E) (33) 

In a circular motion r is the constant and, therefore, the function Af{r) 
has the minimum. Consequently, from the equation dAf/ dr = we find 

/=7V(27 - 3) (34) 

Using (^Uj) we have at f = for E = M : 

e' = 2(7-1)/[7(27-1)] (35) 

From the eqs. ( ^2]) — (|3^) we obtain 

V2^/k/^f3/2l 



Hence, the circular orbits have rotation period 

T = 2V27rc~^a(l + r^f^ (36) 

( 3 rd Kepler law). In comparison with general relativity the change in T is 
2% at r = 3 and 20% at r = 1.33 . 

Let us find the apsidal motion. For ellipsoidal orbits the function Q{u) 
has 3 real roots Hi < II2 < M3 . The apsidal motion per one period is 

5</? = 2|A| -27r, (37) 

where 

A = /_ [g{u)Yl^du (38) 

We have || 

A = 2(M3-Mi)^/2^(7r/2,g), (39) 
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where 

q=[{u2-u,)/{ll,~Ur)Y'^ (40) 

and 

^(7r/2, q) = r'\l - sm\P)Y/'d(3 (41) 
Jo 

Let us introduce (by analogy with general relativity) the following nota- 
tions: 

Ml = (1 - e)/p, 112 = (1 + e/p, % = 1 - 2/p, (42) 

where the parameter p at r — ^ cxd becomes the focal parameter p divided by 
a and parameter e becomes the eccentricity e devided by a . 
At r ^ 1 the value of g ~ (2e/p)^/^ <^ 1 and, therefore, 

^(7r/2, q) = (7r/2)(l + + 9g764 + ...) (43) 

Using eqs. ( ^8l) and (|4ll) we find with accuracy up to l/p^ 

= 3n/p+ {7r/{8f){-e^ + 16e + 54 + ...)• (44) 

For the orbits of Mercury or a binary pulsar (such as PSR 1913 + 16 ) 
the value of u differs very little from the value of a/r . Consequently, the 
values of p = 2/ (ui + U2) and e = (II2 — W/ (^2 + W) differs very little from 
the values of p/a and the e/a . Hence, their apsidal motion differs very little 
from the general relativity prediction. Even, for example, at p = 10 and 
e = 0.5 , the difference in Aip is about 6 ■ 10~'^rad. 

5 Photon Orbits. 

The equations of motion of a photon in the spherical symmetric field are 
given by ||^ 

= (c^ C/A){1 - C b^B), ip = cC b/B, (45) 

where b is the impact parameter. 

The differential equation of orbits can be written as 

du/dip = gi(u), (46) 

where 



13 



and h = h/a. 

Setting r = in egs.® we obtain h = {BCf/"^ = //(I - 1//). Fig 7 
shows 6 as a function of r , i.e. the location of the orbits turning points. 




Fig. 7 The function b = b{r) 



The function b[r) , in contrast to general relativity, is determined on all 
the axes, which is caused by the lack of events horizon. The minimal value of 
b , i.e. bmin ; IS cqual to 2.6 . It is reached at r = 1.5 . The motion of photons 
can be shown by the horizontal b = Const. In Fig. 7 two types of the orbits are 
shown which have the pecularities due to the lack of event horizon. The orbit 
with impact parameter b = bi show that the attracting mass cannot keep a 
photon escaping from the center at the impact parameter b < bmin .This orbit 
type also shows the gravitational capture of a photon. The photon finishes at 
the field center , unlike general relativity, where it end on the Schwarzshild 
sphere. 

The orbit with 6 = 62 shows that the attracting mass can keep a photon 
escaping from the center only if the impact parameter b is greater than bmin 
. The angle of the light deflection at the distances close to r = 1.33 is given 
by 

^ = ln(4.021/52), (47) 
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where 

^ = '^ifmin ~ 3/2) and /^j„ = (1 + T^minY^^ 

6 Conclusion Remarks. 



It follows from eq.(27) that the gravitational energy released by a collapse is 
more than it was supposed to be hitherto. At the same time the gravitational 
potential on the surface of a stable massive configuration of the degenarate 
Fermi - gas is of the order of 

V = {GM/R){1 - exp{-R/a) (48) 

Hence, it follows from the virial theorem that the above objects with R <^ a 
are the ones with low temperatures and, therefore, with low luminosities. 
They,more probably, refer to "dark" matter of the Universe. If their lumi- 
nosities are caused by an accretion, then the Eddington limit of luminosity 
is approximately 

C = C%,[l-exp{-R/a)], (49) 
where ~ 1 ■ lO^^J^Qerg/ s . Hence, if R/a <^ 1 , its luminosity is 

C ^ ^Edd ■ 
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